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We study the nonequilibrium growth of a weakly inter- 
acting homogeneous Bose gas after a quench from a high- 
temperature state to a temperature below the Bose-Einstein 
critical condensation temperature. We quantitatively charac- 
terize the departure from thermal equilibrium and observe the 
presence of two equilibration time scales. The equilibration 
times are shown to be inversely proportional to the density. 
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The nonequilibrium evolution of a Bose-Einstein con- 
densate is a fascinating subject not only on its own but 
also because it allovifs for experimentally probing gen- 
eral ideas concerning the dynamics of phase transitions. 
Since the first experiments with dilute atomic gases led 
to the accomplishment of Bose-Einstein condensation [Q , 
a great theoretical and experimental effort has been put 
in the investigation of different aspects of equilibrium 
and dynamical properties of the Bose-Einstein conden- 
sate To date, however, a comprehensive description 
and complete understanding of the nonequilibrium be- 
havior of the condensate is still lacking. In this context, 
only a few experiments involving the study of the growth 
of the condensate in trapped dilute atomic gases [^[^] 
have been performed. From a theoretical standpoint, dif- 
ferent growth equations, based on the Boltzmann equa- 
tion applied to inhomogeneous trapped Bose gases, were 
proposed ||,^. 

Recently, some of us proposed a method based on the 
closed-time path nonequilibrium field theoretical tech- 
nique to study the growth of the homogeneous Bose con- 
densate 0. We pointed out the fundamental role played 
by the non-condensate fluctuations in the dynamics of 
the condensation growth, and were able to describe the 
condensate as well as the non-condensate components of 
the gas in a unified microscopic formalism without intro- 
ducing any phenomenological parameter. However, from 
a conceptual point of view, the growth of the condensate 
is not the only observable that characterizes the out-of- 
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equilibrium dynamics of a bosonic system. In fact, to 
fully describe the out-of-equilibrium dynamics it would 
be necessary to compute higher-order correlation func- 
tions at two different times. 

Hence, an important issue that deserves a more de- 
tailed study concerns the nonequilibrium process follow- 
ing a rapid quench of the vapor of dilute atomic gases, 
from a state well above the condensation temperature 
to a temperature below critical. Experimentally, this 
can be accomplished by the combined effect of a mag- 
netic trap followed by a rapid RF sweep, as was done 
in the MIT experiment described in Ref. The fi- 

nal equilibrium state of the system depends strongly on 
the nonequilibrium evolution after the quench process. 
Therefore, it is imperative to have a formalism that is 
able to provide not only a description of the formation 
of the condensate but also to quantify the nonequilib- 
rium behavior of the whole system. In another context, 
Gagne and Gleiser [|| have proposed a simple way to reli- 
ably measure the departure from thermal equilibrium of a 
symmetry-broken system which undergoes a temperature 
quench. The method makes use of the absolute value of 
the rate of change of the momentum- integrated structure 
function, AS'tot(i), which we will properly define for our 
case (see Eq. below). The main idea is that for t < 
{i.e., before the quench) the system is in thermal equi- 
librium and AS'tot(i) ~ 0, since in this regime each single 
field mode does not change in time significantly. Simi- 
larly, for large times the system returns to equilibrium 
(with a different temperature) and this quantity again 
approaches zero. Therefore, any departure from zero of 
AS'tot(i) is a measure of the out of equilibrium regime. 
In other words, AS'tot(t) behaves as an "order parame- 
ter" that is sensitive to the out-of-equilibrium evolution 
of the system. In Ref. Q this quantity was evaluated 
for a simple quenched relativistic model and exhibited a 
peaked structure. It was shown that the height of the 
peak measures the departure of the system from thermal 
equilibrium whereas the location of the peak is a precise 
measure of the equilibration time scale, Teq. 

In this work, we apply the formalism developed in Ref. 
to evaluate AS'tot (t) in order to study the nonequilib- 
rium behavior following a quench in a weakly-interacting 
homogeneous dilute Bose gas. The system is assumed to 
be taken from a high temperature equilibrium state, well 
above the critical temperature of Bose-Einstein conden- 
sation, Ti ^ Tc, to an unstable state in a thermal bath 
at a temperature T < T^. Although this may be a strong 
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idealization of the real setup of recent experiments involv- 
ing the evolution of dilute atomic gases subject to a con- 
fining potential, we nevertheless expect that this analysis 
will provide important qualitative insights to the actual 
nonuniform systems. In fact, this approach should gener- 
ically be applicable to the dynamics of trapped atomic 
gases in the central region of wide traps. As we will show 
below, our results reproduce all qualitative aspects of the 
condensate growth. 

Our main result is the determination of a temperature 
range in which there are clearly two distinct equilibration 
time scales. One of the equilibration times is essentially 
determined by thermal fluctuations dominated by pair 
correlations, whereas the other is governed by the con- 
densate growth. We also obtain numerically a scaling 
behavior of these quantities with the density. 

Following Ref. we consider a field theory effective 
model for a nonrelativistic complex Bose field, 0(x, t), 
with mass m, and a hard core interaction potential, 
y(x — x') = g5{x — x'), where the coupling constant, 
(7, is related to the s-wave scattering length, a, as g = 
ATrh^a/m. We use the standard decomposition of the 
fields ((/),(/)*) into a condensate (uniform) part ((/5o,¥'o) 
and a fluctuation (nonuniform) part ((/?, ip* ) that de- 
scribes the atoms (excitations) outside the condensate, as 
<j>{x, t) = (fio (t) + (p{x,t) and (/)* (x, t) = (p'^ (t) + ip* (x, t) , 
where we have assumed a homogeneous condensate. Par- 
ticle conservation enforces the constraint on the total 
density, 

n=\Mt)\^ + {^*^) , (1) 

carried throughout the dynamics. Treating the inter- 
action terms in the self-consistent Hartree-Fock (HF), 
which is also called Hartree-Fock-Bogoliubov approxima- 
tion, we can deduce the following equation of motion for 
the field fluctuations: 

*^"(~^~^) '^ + ^-9"'^ + ^('^oW + (w) (2) 

and an analogous expression for ip* . The quantities {(p(p) 
and {ip*ip*) are usually known as the anomalous density 
terms, while {(p*'p) is the so-called non-condensate den- 
sity. 

In the mean field (or self-consistent HF) approach used 
here, we restrict ourselves to the dynamics in the colli- 
sionless regime (see e.g. Ref. [0), valid when the fre- 
quency of the excitations, Wcxc, is much larger than the 
typical collision (or relaxation) time, icon: '-^oxcicoU S> 1. 
In the opposite regime, Wcxcicoii <C 1, called hydrody- 
namical regime, the dynamics happens in the collision- 
dominated regime in which we can no longer neglect 
three-field contributions like {ip*ipip*) that would enter 
in the equation of motion for the condensate and fluc- 
tuations. These terms describe coUisional effects that 
become important for temperatures close to the critical 



temperature of transition, and must be included in order 
to reproduce the results of the experiments of condensa- 
tion dynamics in a trap. They can be introduced, for ex- 
ample, in the context of quantum Boltzmann equations, 
and lead to a consistent description of the evaporative 
cooling and subsequent formation and equilibration of 
the condensate, leading to a consistent description of the 
experimental data available, as shown in Refs. |^,^. As 
we will discuss below, the parameters we use fall deeply 
into the coUisionlcss regime, warranting the validity of 
the mean-fleld approach. 

In Eq. (|2|), the chemical potential in the HF approxi- 
mation is given by fi^p = g{n+ {f* ip) + {(p(p)) . It how- 
ever does not satisfy the Hugenholtz-Pines (HP) theorem 
for the chemical potential, /ihp — g{n + {^*ip) ~ {w)), 
resulting then in the existence of a gap in the spectrum of 
elementary excitations at equilibrium fl^ , in opposition 
to the known results for Bose condensates. Even though 
we treat the nonequilibrium evolution of the condensate 
where the Hugenholtz-Pines theorem, which is a purely 
static relation valid in equilibrium, is not relevant 
the use of the HF approximation for the description of 
the noncondensate fluctuations violates the conservation 
of particle number, which in equilibrium results in the vi- 
olation of the HP theorem. This is a symptom of the fact 
that this approximation does not respect the constraints 
on the dynamics imposed by the spontaneous breaking 
of a C/(l) symmetry. 

On the other hand, in the Hartree-Fock-Bogoliubov- 
Popov (HFBP) approximation, where the anomalous 
density is neglected in Eq. (|^), the theory is gapless 
p^ . This is the approximation used in Ref. [Q, where 
the chemical potential was shown to he ^ — 2gn. One 
way to go beyond the HFBP approximation and at the 
same time keep the theory gapless, as proposed in Ref. 
jl^, is to treat the (contact) coupling constant g in the 
HFBP approximation as an effective coupling constant 
gcs that can be defined as 

5eff = 5(l + (w>/l'/'on . (3) 

As shown by the authors in Ref. jlj] , this is a consistent 
way to include the effects of pair correlations, represented 
by the anomalous density, and accounts for the low mo- 
mentum limit corrections of the many body T-matrix 
keeping the theory gapless. 

Eq. (||) in the improved HFBP approximation together 
with the constraint on the total density completely de- 
termines the time evolution of the condensate, the non- 
condensate density {ip*ip), as well as the anomalous den- 
sity (ipip) in ^off . In order to solve these equations numer- 
ically, it is convenient to rewrite them in terms of their 
real field components, ip ~ tpi +iip2 and its complex con- 
jugate. The two-point function {ip*'p>) and the anoma- 
lous densities can be expressed in a convenient way in 
terms of the Green's functions (tpiipj) {i,j = 1,2) con- 
structed from the homogeneous solutions to the operator 
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of quadratic fluctuations 
Keldysh formalism (see, e.g. 



Using the usual Schwinger- 
find § 
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{LPj{t)Lpj{t)) = 



(2^)3 1 - exp (-/3£k) 
+ exp(-/3ek)x*(k,t)x,(k,t')] 



[x,(k,t)x*(k,t') 



(4) 



where £k = k^/(2m), and (3^^ is the final temperature 
the system will reach at long times. The equations for the 
fiuctuation modes, xi and X2, that are the homogeneous 
solutions to the operator of quadratic fluctuations, were 
derived before in the HFBP approximation. These 
mode equations, in the improved HFBP approximation 
discussed above, become 
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XI = 
X2 =0 
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where g^s is given by Eq. (^). The boundary con- 
ditions for the solution of Eq. (||) are such that for 
t < 0, |iy5o(OP = and XiiKt) and X2(k, are 
given by xi(k, i) = exp(iekt)/\/2iV and X2(k, i) = 
i exp(ieki) /V'2N . The normalization factor N is easily 
found to be § iV = 2{Po/Pf^'^ where is the critical 
temperature for condensation of an ideal Bose gas. 

The density constraint, when expressed as a solution 
for the condensate density no{t) = \ipo{t)\'^, becomes 
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|X2(k,i)H]nk(/3) , (6) 



where nk(/3) = {e^^^ — 1)^^. The momentum integral 
in Eq. (|^) is bound to values of fc < /cc, corresponding 
to the unstable modes that drive the system from the 
initial excited state to the condensate state, k^ can be 
determined from the equation for the modes, and is given 
by the value of k which changes the sign in the product 
of the frequencies in Eq. (|^). 

In Fig. 1 we display the result for the evolution of the 
condensate density, as given by the solution for the cou- 
pled integro-differential system of equations (||) and (|^), 
for three values of the final temperature. This result im- 
proves over the previous one, presented in Ref. j^], which 
was obtained in the HFBP approximation by neglecting 
the contribution from the anomalous densities. Although 
the HFBP approximation describes well the condensate 
for temperatures well below the critical temperature, it 
is well known that it breaks down for temperatures of 
the thermal bath near the critical one, in which case the 
anomalous densities become of the same magnitude as 
the out of the condensate excitation densities. In partic- 
ular, we show in Fig. 2 that even for moderate temper- 
atures the anomalous densities can grow to appreciable 



values during the out-of-equilibrium evolution of the sys- 
tem. 

From the results of Fig. 1 we have verified that, at large 
the temperature dependence of the condensate (in equi- 
librium) is the expected one, matching the known results 
for the drop of the condensate density in equilibrium for 
increasing temperatures in an interacting gas How- 
ever, for initial times, this tendency changes due to pair 
correlation contributions represented by the anomalous 
density, in Eq. (||). This fact suggests the exis- 

tence of two different scales for the equilibration time: 
one related to the condensate growth and the other to 
the anomalous density. 

In order to investigate this behavior more precisely, let 
us consider the structure function defined as the Fourier 
transform of the two-point Green function, S'fc(i) = 
{(pj{t)(pj{t))k [0 (this should not be confused with 
the usual dynamical structure function, given in terms 
of the Fourier transform of the four-point Green func- 
tion instead). We find it more convenient to work 
with the normalized structure function, Sk (t) / Sk (0) , 
which, using Eq. (||), is given by Sk{t)/Sk{0) = 
iV (|xi(k,t)p |x2(k, i)P). The wave-number inte- 
grated time derivative of the normalized structure func- 
tion provides the net change in the excitations outside 
the condensate, or the total change of fiuctuations in the 
condensate. We define this quantity as 



A^tot(0 = 



^'4 



Skjt) 
Sk{0) 



(7) 



As shown in Ref. g, the peak of AS'tot(0 gives a mea- 
sure of the departure from thermal equilibrium following 
a quench, while the position of the peak gives the equi- 
libration time scale of the system. These properties can 
be extracted from Fig. 3, which shows AStot{t) for the 
three cases defined in Fig. 1. 

In Fig. 3 we observe two peaks for each temperature 
studied. The first peak is clearly related to the time evo- 
lution of the anomalous density (compare with Fig. 2), 
and the second peak to the time evolution of the conden- 
sate density (compare with Fig. 1). Notice the correla- 
tion of the results regarding the time scale of equilibra- 
tion for the condensate (associated to the position of the 
second peak). The lower the temperature of the ther- 
mal bath into which the system is quenched the higher 
is the peak of AS'tot(i)- This is consistent with the fact 
that the lower the temperature of the thermal bath the 
further from equilibrium the system will be. The tem- 
perature dependence of the equilibration times is quite 
clear. For low temperatures, T/Tq < 0.1, the first peak 
disappears. The equilibration process in this regime is 
completely driven by the condensation growth. In this 
case, the approximations used in Ref. j^] are acceptable 
since the anomalous density contribution becomes sub- 
leading. On the other hand, for final temperatures near 
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the critical one, 0.3 < T/Tq < 1, the second peak in Fig. 
^ disappears and the equihbration process becomes now 
fuUy governed by thermal fluctuations encoded in the 
mean value of the anomalous density. Moreover, we find 
an intermediate temperature region, 0.1 < T/Tq < 0.3 
(shown in Fig. ^), where the two equilibration times are 
present. 

In Fig. 4 we exhibit the results for the equilibra- 
tion time tcq (= {ma^/h)Tcq) for both the condensate 
(squares) and the anomalous density (triangles) as a func- 
tion of the (dimensionless) density p = na^ at fixed tem- 
perature T/Tq — 0.1, extracted from the location of the 
second and first peaks, respectively, of AS'tot. The results 
can be perfectly fitted by the function t^^^2 — 0.2/p 
(line), for the condensate growth, while the result re- 
lated to the equilibration of the anomalous density scales 
with the density in an analogous way, t^^ ^^^s^ — 0.12/p 
(dashed line). In terms of dimensionful quantities: 



These results could in principle be compared with some 
early qualitative estimates for the growth time of the 
condensate for homogeneous gases [0 , even though the 
physics discussed in those references is essentially dif- 
ferent, yielding a rough estimate of the formation time 
~ 0[{ng)~^]. This is exactly the behavior found here. 

Although we deal with an idealized situation for the 
condensation formation, it is useful to briefiy discuss the 
experimental results and compare them with our esti- 
mates, pointing out clearly the differences and similar- 
ities. Up to date only two experiments have been per- 
formed involving the explicit study of the dynamics of 
BEC formation. The first one was performed by the 
MIT group and reported in Ref. H]. They have stud- 
ied the dynamics of BEC of atomic ^^iVa (a — 2.75 nm) 
in a cigar shaped magnetic trap, described by frequen- 
cies ujx — ojy ^ 27r82.3 Hz and ojz — 27rl8 Hz. The 
number of atoms in the trap was N ~ 10^ and the 
temperatures involved were between 0.5 — 1.5 fxK. The 
geometry of the confining potential yields a volume at 
the center of the trap ||^ Vtrap = Tr^^^axOyaz, where 
Ui — Ti / (muji) . The central density can then be es- 
timated as n ^ 6.8 x IQ-^^ atoms/m^, corresponding to 
p = na^ ~ 10"'^. The typical formation time observed 
was - 100 - 200 ms. 

In the second experiment, Ref. performed by 
the Munich group, they used (a — 5.77 nm) in 

a trap characterized by ujx = c^y — 27rllO Hz and 
LUz — 27r 14 Hz. The number of atoms in the trap 
was N ~ 4.2 X 10^ and the temperatures involved was 
~ GAOnK. These parameters lead to a central density 
value of rt ^ 2.5 x 10^'^ atoms/m"^, corresponding to 
p = na^ ~ 5 X 10^^. The typical formation time ob- 
served was also ~ 200 ms. 



From the experimental data above, we can estimate in 
which regime is each of the experiments. For trapped 
atom gas experiments Woxc scales with the trapping po- 
tential frequency and tcoi is usually a measured quantity. 
Taking (jJoxc ~ w = {uj^tOyLOzY^^ , the average frequency, 
and using the experimental data for the relaxation time 
measured in these experiments, we find for the MIT ex- 
periment that Wcxc^coU ~ 3, while for the Munich exper- 
iment ti-'cxc^coU ^ 40, putting these experiments in the 
borderline between the coUisionless and hydrodynamical 
regimes (MIT) and in the coUisionless regime (Munich). 
For a homogeneous gas it is more difficult to get such esti- 
mates. However, it is reasonable to estimate cjexc as pro- 
portional to the thermal energy of the bath ?iWexc ~ ksT 
and tcoii expressed in terms of the mean time between 
collisions, which for a classical gas is {anv)~^ (for an ex- 
plicitly determination of the relaxation time, as given by 
the inverse of the damping rate for a homogeneous weakly 
interacting Bose gas, see Ref. [0), where a = Stto^ is the 
s-wave cross-section, n is the maximum non-condensate 
density and v is the characteristic thermal velocity. Using 
our parameter values and a temperature of T = O.lTc, 
with Tc approximated by the condensation temperature 
for an ideal gasQ, we obtain Woxc^coii ~ 60, deeply into 
the coUisionless regime, for both '^^Na and ^"^ Rh. The 
results for the equilibration times we obtain are of order 
lO^^s for the parameters of ^"^ Rb and lO^^s for "^^Na. 

The large difference in magnitude of time scales ob- 
tained here as compared to those of the experiments sig- 
nals the importance of incoherent /coherent coUisional ef- 
fects in the trapped atomic gases experiments due to the 
nature of the experimental setup. Differently from our 
quenching scenario, where we can effectively approximate 
the model by a two-level system at low temperatures, in 
the experiments, the confined gas is firstly cooled to a 
temperature close to the critical temperature and then it 
is put in a nonequilibrium state by means of an evapora- 
tive cooling through a rapid rf sweep (MIT experiment) 
or a constant frequency rf field (Munich experiment), 
which removes the highly excited atoms from the trap. 
The system in then left to thermalize to a temperature 
below Tc and the formation (growth) of the condensate 
occurs. This cannot be approximated by a two-level sys- 
tem. In the experiments, scattering and growth among 
adjacent levels, which are described by those microscopic 
processes associated with collisions, are of relevance to 
the dynamics in this regime. Those processes are in- 
deed appropriately described through the use of quan- 
tum kinetic theory (quantum Boltzmann equations) for 



For an estimate this is a reasonable approximation since for 
a dilute weakly interacting gas the the shift of the transition 
temperature, ATc/Tc, of the interacting model compared with 
the ideal gas transition temperature, is of order an}!^ [^. 
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the trapped gases and are necessary to describe the whole 
process of evaporative cooling and subsequent thermal- 
ization during the process of condensate formation as 
shown in Refs. It is clear that we are considering a 

different regime from the beginning. In particular, we do 
not have a thermalization regime (since it is considered in 
thermal equilibrium from the beginning). CoUisional ef- 
fects lead to a much longer evolution for the condensate, 
explaining the disparities of orders of magnitude between 
our results and the experimental ones. We recall that, as 
pointed out by Stoof in Ref. jl3j, the regime we stud- 
ied (named the second stage of condensate formation in 
that reference) lasts about 0{h/ng), which corresponds 
exactly to the time scale we found, while the relaxation 
towards equilibrium in the third and final stage of the 
condensate formation (absent in our theoretical setting) 
leads to a much longer time evolution for the condensate. 

In summary, we have provided in this paper a way to 
analyze the time scales related with the formation of the 
condensate, by studying the behavior of the total struc- 
ture function, as defined by Eq. (Q). We have analyzed 
the case of a homogeneous gas tha undergoes an instan- 
taneous quench and we have shown that for tempera- 
tures of the thermal bath much lower than the critical 
temperature the gas is deep in the coUisionless regime, 
justifying the mean field approach developed in this pa- 
per. We showed that under this circumstances, there are 
two equilibration time scales perfectly distinguished and 
associated with different physical processes, i. e. the 
growth of the condensate and the pair scattering process 
described by the anomalous density. The structure of 
equilibration times reported in this paper opens several 
questions related to the approach to equilibrium in in- 
teracting dilute gases. For instance, the presence of two 
equilibration time scales may infiuence time correlations, 
modifying fluctuation-dissipation relations. In order to 
address these issues, theoretically as well as experimen- 
tally, it is necessary to understand the behavior of corre- 
lation functions that depend on two different times as, 
for instance, the two-time dynamical structure factor. 
We hope to report on this subject in the near future. 
Another issue we are currently pursuing is the improve- 
ment of the equations presented here by including the 
higher-order collision terms by means of self-consistent 
Schwinger-Dyson equations for the nonequilibrium prop- 
agators. These results will be presented in the near fu- 
ture. 
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FIG. 1. Fraction of condensate density as a function of time 
for no? = 0.01 and T/Tq = 0.1,0.2,0.3. Here, t = {h/ma^)t 
is a dimensionless time. 
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FIG. 2. Time evolution of the anomalous density for the 
same parameters used in Fig. 1. 
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FIG. 3. Time evolution of the wave number integrated time 
derivative of the normalized structure function for the same 
parameters used in Fig. 1. 
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and anomalous density (triangles) obtained from the second 
and first peaks of AS'tot, respectively, as a function of the 
density, p = na^ , at a fixed temperature T/Tq = 0.1. 
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